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(N , 

CN ! This is the second of the two last papers by V. N. Gribov concluding his 20 year 

long study of the problem of quark confinement in QCD. In this paper the analytic 

Q^ ■ structure of quark and gluon Green's functions is investigated in the framework of the 

theory of confinement based on the phenomenon of supercritical binding of light quarks. 

Q^l The problem of unitarity in a confining theory is discussed. The write-up remained 

o ■ unfinished and so it is presented here. The author was planning to emphasise the link 

O , between the electro-weak and strong interactions, and in particular the role of pions 

(Goldstone bosons) in confinement, to present an explicit solution for bound states, and 

to write down an analytic model for quark and gluon Green's functions corresponding 

^ ■ to confinement. 



1 Introduction f] 

Almost ten years ago [1] I proposed a hypothesis according to which quarks and gluons are 
confined due to the existence of light quarks {u and d) with Compton wavelengths much 
larger than the radius of strong interaction defined by Aqcd- The essence of this hypothesis 
is the supercritical phenomenon which is well-known in QED where it is the following. If 
there exists a heavy nucleus of radius R with charge Z exceeding a critical value Zcr (which 
is of the order of 137), the vacuum of the light charged fermions (electrons) becomes unstable 

* The work was completed in Bonn under the Humboldt Research Award Program. 

^ The text was prepared for publication by Yu. Dokshitzer, C. Ewerz, and J. Nyiri, with help from 
J. Bartels, A. Kaidalov, A. Mueller and A. Vainshtein. 

-'^The author intended to rewrite this Introduction. In the present form it corresponds to the Introduction 
to the Lecture given by V. N. Gribov at the 34*'* International School of Subnuclear Physics in Erice, Italy, 
in 1996 [7]. 



due to the process of (e~, e+) pair creation. The negative component of the pair, the electron, 
falls into the heavy charge and the positive component goes to infinity. The condition for 
this phenomenon to occur is that the Compton wave length 1/m of the electron has to be 
much larger than the radius of the heavy charge: 



a 



v/^2 _ zl ln(mi?) > 71 . (1.1^ 



The electron which falls into the centre forms together with the source Z an ion of charge 
Zi = Z — 1. li Z — 1 > Zcr, this process will continue until the vacuum becomes stable i.e. 
when Zn = Z — n becomes less than the critical charge. 

From the point of view of the Dirac equation the falling electron has negative energy. On 
the other hand, the supercritical ion with Zn < Z^ is stable because of the Pauli principle. 
Indeed, the electron cannot leave the ion, since all negative energy states outside it are 
occupied. 

In QCD this phenomenon can occur since, due to asymptotic freedom, the colour charge 
can reach the critical value for any object. Let us suppose that two heavy quarks are created 
with opposite colours. The gluonic vacuum polarisation increases the colour charges of the 




Figure 1: 

quarks (see Fig. 1). When these colour charges become large enough, light antiquarks start to 
fall onto the heavy quark. However, the bound state formed by a light antiquark and a heavy 
quark will be very different from the bound states which appeared in QED. In QED we had 
an ion with a charge Z — \. In QCD the corresponding state can be colourless. The difference 
is due to the fact, that while the large colour charge of the heavy quark results from gluonic 
vacuum polarisation, the bound state is formed by two particles the total charge of which 
equals zero, and this charge cannot be changed by vacuum polarisation. This phenomenon 
can be called quantum screening. Formally it means that the heavy quark q^ is unstable and 
has to decay into a meson and a light quark qi 

qh^M + qi. (1.2) 

This case — the confinement of the heavy quark — is relatively simple and almost inde- 
pendent of the properties of the light quark [2]. The problem of the light quark is more 
difficult. At first sight it is not clear whether the "falling into the centre" and the formation 
of a supercritical state can occur in the system of a light quark and a light antiquark. It 
was shown [3], however, that indeed the critical phenomenon exists in the light qq system. 



Moreover, the critical coupling constant is proven to be sufficiently small: 
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There is an even more complicated question: what kind of states can be formed as a result 
of this phenomenon? If the appearing new state is a normal meson, the "falling into the 
centre" leads to 

qi^M^qi (1.4) 

instead of (|1.2| ). Equation (|1.4| ) can be satisfied if there are not only positive energy quarks 
g_i_ but also negative energy quarks g_. This result, although it looks very strange, is not 
unexpected. Indeed, we already learned in QED that particles which fall into the centre 
correspond to the negative energy solution of the Dirac equation. 

To resolve the problem of the negative energy states, Dirac supposed that they are all 
occupied and therefore not observable. The necessity to consider now both the positive and 
negative solutions of the Dirac equation means only that Dirac's hypothesis is not always 
true. If the interaction is strong enough, there is another possibility: negative energy states 
might be only partly occupied and positive energy states only partly empty. This is, however, 
not in contradiction with the absence of negative energy particles in the real world if stable 
supercritical bound states (mesons) exist since in this case both the positive and negative 
energy quarks will be unstable: 

g+ ^ M + g_ 

q_^M + q_q_q_ (1.5) 

and, consequently, unobservable. 

From the point of view of the fermionic spectrum in the vacuum, ( |1.5|) means the follow- 
ing. In the case of weak interactions, the fermionic spectrum has a structure corresponding 
to Fig. 0. All levels above go = \/q^~+~rn? are empty, and the levels below go = — A/g^ + m^ 
are occupied. If the interaction is stronger than critical, quark-antiquark pairs are present in 
the vacuum, and the spectrum has a structure which roughly corresponds to Fig. ^: all levels 
are occupied in the shaded regions and empty in the unshaded ones. In terms of condensed 



matter physics, the curves go = itwg*^ + m^ correspond to the Fermi surface. The change in 

the number of constituents in the theory is also natural from the point of view of the physics 
of condensed matter. In a theory with weak coupling we had two constituents g and g; a 
theory with supercritical coupling contains four states g+, g_, g+ and g_. In non-relativistic 
physics there are only particles (electrons) and no antiparticles. But in a conductor we have 
particles and holes, and holes are negative energy states (counting the energy of the holes 
from zero and not from the Fermi energy as it is usually done). 

In a relativistic theory it is impossible to add a constant to the energy of particles, and 
therefore we are forced to talk about negative energy states which exist only inside our 
matter (vacuum). 
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Figure 3: 



This unusual quark spectrum has been discussed for several years. Still, I have not been 
able to formulate a constructive theory because I have not understood what meson had to 
be introduced as a supercritical bound state. From the structure of the spectrum shown 
in Fig. 1^ it is clear that there have to be different types of excitations corresponding to 
different types of mesons. It is natural to identify hole-particle excitations of the type 1 
near the Fermi surface with mesons like p, u etc. Hole-particle excitations near the light 
cone can be identified with the Oq and /o mesons [4] and also with rj'. But all these mesons 
do not look like natural candidates for the lowest supercritical bound states. Only recently 
did I recognise that there are strong reasons to believe that the pseudoscalar octet, and the 
TT-meson in particular, are in fact the lowest supercritical bound states. Most of these reasons 
are connected with the discussion of the nature of the vr- meson state which I presented in [5] . 

As I have said before, the supercritical atom in QED is not a usual Bohr- type atom, 
containing a definite number of electrons. Rather, it is a collective state in which only the 
electron density is localised around the heavy charge, and only at large distances (much 



larger than the atomic radius) it looks like a state with a definite electron number. A quasi- 
Goldstone state like a vr-meson with a finite mass has the same properties. At short distances 
(distances less than Aqcd) it can be considered as a collective state of the type qiqi — qrQr 
(divergence of the pseudovector current; qi and g^ stand for left-handed and right-handed 
quarks), whereas at large distances it is the two-particle state qiqr — qrqi- Accepting this 
identification with the vr-meson as the lowest supercritical state by introducing it explicitly 
in the equation for the quark Green function (which will be discussed in the next part of the 
lecture), it proves possible to find a solution for this Green function which has properties 
corresponding to the spectrum in Fig. ^ This solution has two complex poles as functions of 
energy in the complex energy plane in accordance with the two types of quarks (constituent 
and current) we have discussed. It has no pole on the real axis which guarantees that quarks 
as propagating states do not exist. It has a soft singularity when g^ — > 0, refiecting the 
fact that quark currents exist in the region where the hadrons are created. It is important 
to stress that in all these considerations I assume the coupling constant a is saturated at 
a value not much larger than acr- Quark masses corresponding to the positions of the two 
complex poles m± are of the order of 

m^ ~ Aqcd exp 



^/a — ac 

The solution for the Green's function is self-consistent if the vr-meson mass is close to m_. 
This means that at least in the case when a is not very large, the vr-meson mass is defined 
by strong interaction dynamics. This result contradicts the usual point of view according to 
which the vr-meson mass squared is the product of the bare quark mass rriQ defined by weak 
interactions and the condensate density {^^) / fl defined by strong interactions. 

In what follows I will explain what type of equations have to be written and solved in 
order to come to the conclusion I have stated and to calculate quantities which have not 
been analysed until now. For example, quark-gluon vertices and the gluon Green's function 
have not been calculated yet. It is clear, however, that the gluon Green's function we obtain 
will also have a complex singularity due to the gluon decay into a qq pair. Consequently, 
the gluon will also be confined. Equations for hadronic amplitudes can also be written 
constructively if a is not very large. It can be shown that if the quark Green's function has 
properties as described above, the hadronic amplitude will have no singularities connected 
with intermediate qq states, but will have singularities related to the vr-meson thresholds. 
The fact that the vr-meson is included intrinsically in the equation for the quark Green's 
function and has a mass of the order of m_ makes this statement much less mysterious than 
it would look without it. 



2 The structure of the confined solution for the 
Green's function of massless quarks 

In the paper [6] we considered the solution for the quark Green's function corresponding to 
the chiral symmetry breaking. We discussed the importance of including the contribution 
of the Goldstone boson in the equation. In the present paper we show that this solution 



does not necessarily survive in the presence of the Goldstone contribution which, essentially, 
reflects the softness of the condensate and leads to the existence of a solution which has no 
poles and which corresponds to confined quarks. 

We will accept the equation for the quark Green's function in the form 



d'G'\q) = g{q)d,G-\q)G{q)d,G-\q) - 
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{^l^.G-\q)}G{^^,,G-\q)} (2.1) 



where g{q) = |^ is supposed to behave like in [6] (see Fig. ^ and /^ is the amphtude for 




Figure 4: 



the pion - axial current transition. 
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(2.2) 



Before turning to the formal solution, let us discuss its general properties and the difference 
in the structure caused by the inclusion of the pion contribution. 



Writing G ^{q) as 



G-\q) = Z-\q){m{q)-q), 



(2.3) 



we obtain (|2.1| ) in the form 



d^G-' = g{q)d^G-'Gd^G-' + 



m m — q 

J2 ^ _ g2 



z-\ 



(2.4) 



where /^ = |7r^/^. For the solution of the equation for G^^ without pion contribution (i.e. 
without the last term in ( p.4|) ) the behaviour of Z{q^) and m{q^) as functions of q^ in the 
Euclidean region of negative q^ = —Q^ is shown in Fig. ^. 

Here the dashed curve corresponds to the massive quark, ttiq 7^ 0, and the solid curve to the 
massless quark, respectively. In the latter case 
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Figure 5: 

If we include the last term of (|2.4| ), the solution corresponding to massive quarks disappears, 
as it was expected. The massless quark solution ( |2.5| ) will, however, survive because in this 
case the last term in (|2.4| ) which can be called A^r is small at large Q^ values: 
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(2.6) 



Let us see what would be the effect of the pion contribution at small Q^. For Q <^ uic = m,{0) 

(2.7) 



4]) leads to the following two equations: 

d^Z-'m) = 
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(2.8) 



We see that the pion produces a hundred percent correction to the effective coupling because 
/ ~ m,c and g ~ 0.2 if we are close to the critical value of g. This change can be crucial. 
We shall analyse the equation ( p.6| ) more carefully later. For the time being, in order to see 
what can happen, let us accept (|2.71 ) and (|2.8|) literally. These simple equations are easy to 
solve. The solutions are 



(2.9) 
(2.10) 



where Zi and Z2 are solutions Z^, of the Bessel equation with index z/ = 1, 2, respectively. The 
concrete forms of Zi and Z2 depend on the boundary conditions imposed on the solutions. 
If we want to preserve the behaviour of Z~^ and m at small Q^ values, we have to choose 
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(2.11) 



with J\ and J2 the Bessel functions. In this case, however, the behaviour oi Z ^ and m, will 
change at large Q"^ values and it is not obvious whether it will be possible to preserve the 



asymptotic behaviour (Fig. |^) corresponding to asymptotic freedom. If our aim is to keep 
the behaviour at large Q^ unchanged we have to choose Zi and Z2 as superpositions of the 
singular and non-singular solutions of the Bessel equation, 

aiFi(Q//ii) + 61 Ji(g//ii) , (2.12a) 
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a2>"2(Q//X2)+&2^2(Q//i2) 



(2.12b) 



and to select 01^2 and hi2 such that the asymptotic behaviour at large Q^ is preserved. For 
small Q^ values this means 
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(2.13) 



(2.14) 



Instead of the behaviour oi Z ^ and m corresponding to Fig. ^ we now have a behaviour as 
shown in Fig. ||. 




Figure 6: 



This behaviour corresponds to the confined solution for the quark Green's function. In this 
solution the condensate which is created at momenta of the order of A exists only in a region 
of Q^ values between Ql and Q| {Ql, Q2 ~ /^) ^^^ disappears at smaller Q^ values due to its 
decay into vr-mesons. Because of the decrease oim{q'^) at large and small Q^, in this solution 
the pion contribution to the right-hand side of the equation (P^), A^, is localised between Ql 



and Q2 (we shall see this in detail later). This localisation of the pion contribution enables 
us to analyse not only the behaviour in the Euclidean region g^ < but also the analytic 
properties of the solution and its behaviour at positive g^ values. 

If the pion is localised, G^^{q) satisfies the old equation without A^ at Q^ <^ Ql and 
Q^ ^ Q2- This means that G^^{q) can be written in the form 

G-\q) = CMG^\q) + C2{q)G2\q) , (2.15) 

where G^\{q) are two different solutions of the equation ( p.4| ) for A^r = and 6*1^2(5') are 
slowly varying functions. 



Ci{q) -^ const at Q^ :^ Ql and Ci(g) ^ at Q^ <^ Q 
G2{q) -^ at Q^ :$> Ql and ^2(5) -* const at Q^ <^ Q 



2 
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G~^{q) has to satisfy the conditions at large Q"^ which we already imposed on our solution. 
For this, we are bound to choose Gi^{q) to be the solution corresponding to symmetry 
breaking with the "mass" rric, which was discussed in [6], and G^^(g) has to be the solution 
describing a singular behaviour of the type ( p.l3| ), ( p.l4| ) at small Q^. We have to have in 
mind also that G^^{q) (|2.15|) must not have singularities in the complex plane. G'j~^(g) has 
no singularities there but the standard cut on the real axis from g^ = ml to g^ ^ oo. We 
choose G^^ to have also no singularities in the complex plane. After that we shall see what 
will be the singularities of G~^{q). 

Making use of the fact that equation (2^) with A^ = is scale invariant for a = const, 
we can always write 
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(2,16) 



where G ^(g) satisfies the equation 



m\ 



d'G-\q) = W Y ) d,G-\q)G{q)d,G-\q) . 



(2.17) 



This is the same equation which we have discussed in [6]; the behaviour of g{q) = girn^/q) 
is ^ ^ at g ^ and ^(g) — > (?o at g — i> oo, opposite to that shown in Fig. §. It was shown 
in [3] and [6] that independently of the behaviour of g we can always find a solution with a 
cut along the positive g^ axis from g^ = m^ to g^ — >■ oo if we choose the solution not to have 



singularities at g = and if we fix G ^(g) by the condition G ^{q)\q= 
in a simple relation). 
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Suppose that G ^ (g) is chosen in such a way. Then G2 ^ (g) defined by (|2.16|) will have a 



cut from g = to g = mt 
cuts in the g^ plane (Fig. |^), 



and a singular behaviour at g^ ^0. As a result, G ^(g) has two 
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Figure 7: 



where we take m,2 > m,i since, as we will see in the next section, this is the only possibility 
to avoid singularities in the complex plane; due to reasons to be discussed, ie is positioned 
in the way shown in Fig. 0. 



For the sake of simplicity, let us consider the properties of G ^(g) in the go-plane at the 
value g = of the space component of q. We will write G~^{q), as in [6], in the form 



G-\q,) = G7^(go)^ + G-_\qo)^-^ 



(2.18) 



For example, ^+"^(^0) has two normal cuts (Fig. 
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Figure 8: 



The singularities of Gi+(go) at go = "^i and go = 
(103) in [6] that Gr+(9o = rrii) = 0, while ^^^(go 
Gi^iqo = mi) = Gil{qo = —mi) = const. 



-mi are different. It is clear from eq. 
mi), for (7 < |, is different from zero: 



The function G'^^^(go) also has the analytic structure of Fig. ^ with mi — ^ m2- In the 
course of the reflection go — > —rn^/qQ the points on the line go = iQ transform into gg = 
im\/Q] the points go = ±m2 =F i^ transform into go = =F^2 =F i^i respectively. As it follows 
from all this, G^^ has singularities corresponding to Fig. Rl. 
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Figure 9: 



These analytic properties are a clear manifestation of the fact that the Dirac see is 
destroyed; G^^ has singularities corresponding to both positive and negative energies. The 
parameter m2 has the meaning of the Fermi energy. In order to find zeros of G^ which have 
to be in the lower go half plane and correspond to unstable quarks with positive and negative 
energy, we have to know the signs of ImGf.i';, IniG^.^ and the signs of Ci, C2 in (|2.15 ). The 
solutions Gi"_^(go) and G^{qo) satisfy the normal unitarity condition. Their imaginary parts 
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at positive go have to be negative. This leads to signs of the imaginary parts as shown in 
Fig. 1^ (with arrows pointing at the positive side of the cuts). Having this in mind, we can 
write Gi^{q) near to its zero (i.e. at qo close to mi) in the form 



G^l = {nil -iS -go)" 



<(3<1. 



(2.19) 



Since 6*2+ remains finite at go -^ mi, it can be given in the form 



G^Uqo 
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(2.20) 



If so, 



G';^(go) = Ci{mi -i5- go)^ + Cape-*^. 



(2.21) 



The zeros of Gj^ (go) are then defined by the equation 



(2.22) 



If ^ > 0, we have ImgQ < 0. In this case the singularity appears in the lower half plane 
and describes the unstable positive energy quark. Repeating the calculation for go — > — m2, 
we obtain a singularity corresponding to the unstable hole. Considering Gl^(go) instead of 
G^{qQ), we find the same singularities for an anti-quark and an anti-hole. 

If the condition ^ > is not satisfied the singularities can move to the upper half plane. 
This does not destroy the theory because they will be on the unphysical sheet. However, in 
this case I do not have any simple interpretations for these singularities. 

The analytic properties of G^{q) presented in Fig. ^ imply that in the limit e — ;> the 
Feynman Green's function defined for real go values becomes a non-analytic function. Two 
analytic functions which have no singularities in the upper and lower complex half planes 
can be defined and called the retarded and the advanced Green's function. At finite ie there 
exists one analytic function with four cuts. This can be easily seen if we consider G^^{q) as 
a function of go at a fixed value of the space component g*of the 4- vector g^. In this case we 
will have Fig. |10| instead of Fig. ||. 
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When negative energies are involved, it is natural to expect that for the exact solution all 
four cuts have discontinuities different from zero in the intervals from go = 1^1 to go ~^ cx> 
and from go —>■ — oo to go = — |g |. 

If in the usual, non-confined case we know the Green's function (calculating it by using 
Feynman diagrams in the Euclidean space), we can continue it into positive g^ values and 
find that the discontinuity at g^ > m^ satisfies the unitarity condition. This means that the 
retarded Green's function coincides with the Feynman Green's function. 

The confined case is more complicated. The knowledge of the retarded Green's function 
in the limit {ie — > 0) is not sufficient for finding the Feynman Green's function. In order to 
obtain the Feynman Green's function, we need the equation for the discontinuities on the 
new cuts; it will be an equation for the density matrix of quarks in the vacuum. 



In the next section we will find the solution of the equation (|2.1|) in the limit ie -^ 
(retarded Green's function) without poles and with the properties we have discussed. In 
section ^ we shall obtain the equations for the discontinuities on the new and old cuts. 



3 Solution for the retarded Green's function 
of confined massless quarks 

In this section I follow the pattern formulated in the previous section. I will find the solution 
of the equation ( p.l|) in Euclidean space and then, to be sure that it is stable, continue it 
into the complex plane. To find the confined solution of the equation ( |2.1| ) for the Green's 
function of light quarks, we will introduce the same representation for G~^ as in [6]: 

G-fe)=(^)%--l (3.1) 

where (3 = 1 — g and h = q/q. Instead of the equations (84), (85) in [6], we will have 

2u ■ 
0+— 0-3sinh0 = O, (3.2a) 

u 



u — u + 



,M3si,>h=^4 -?«i 



u = —{ii — u) , (3.2b) 

P 



where / = gi and ^ = Ing. As we see, the pion contribution to the equation (O) influences 
only the equation for m(^) and it depends explicitly on g^. Because of this, for the equations 
( p.2|) the energy is not conserved even if /3 is constant. In the Euclidean region, g^ = —Q^ < 0, 
we can write (p = iip and the equations take the form 

t/^ + 2pip — 3 sin-?/; = , (3.3a) 

u - V{^)u = , (3.3b) 
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where 



V'K) = l+,a^|3sin^| + f)-^floos^| + i(p 



(3.4) 



and p = ii/u. As before, the equation (p.3a|) corresponds to particle propagation in a periodic 
potential with damping when ^ is increasing and with acceleration when ^ is decreasing. 
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Figure 11: 



In [6] we considered the trajectories of the type I and II (dashed curves); now we will 
concentrate mainly on the trajectory C (solid curve) which has the structure we discussed in 
the previous section. Any given trajectory defines the potential V{^) (if we neglect ^(ii — u)) 
in the equation (|3.3b| ) which is the Schrodinger equation at zero energy. The solution of this 
equation, defining the damping p in ( p.3a| ), has to be chosen self-consistently. 



The structure of the equations ( p.3|) with pion contribution differs essentially from the 
structure of the equations without pions even at large Q^ in the sense that they have no 
solutions corresponding to massive quarks. The massless solutions of these equations are not 
very different with or without vr-mesons. 



The solutions of 



-(ip-n) 
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satisfying the condition of asymptotic freedom are, at large Q^ 
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where 7 = | is the invariant anomalous dimension; b = yA^c — |^/- The renormalisation of 
the Green's function which we discussed in the previous section is 



7//3 



Z-\Q)=uo{ — ) at Q 

in Feynman gauge (which we are using here). 
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(3.6) 
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It is not clear at all whether it makes sense to take these anomalous dimensions seriously. 
For the sake of simplicity, we shall consider ( |3.5| ) at 7 = 0. Important is, that the solutions 
of ( |3.3|) still contain two parameters. 



For small Q^ values the non-confined solutions can be written in the simple form 

Q 



ij 
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m. 



Q{i + lm(m$^~^ 



The confined solution has a more complicated structure. For Q^ 
this case 
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(3.7a) 
(3.7b) 
0, ip ^ TT we have in 

(3.8a) 



u 



+ f-^V'-££L cos (2^/2^:3:5^ In ^ + (5)1, (3.8b) 
\m2j ^2 - 3/32 ^ "^ "^ go J 



with p = ^/l+W, P = /3(0). 

The solutions contain three essential parameters C, Qq and Qi which define ip = ip — ir , ip 

and -it at Q = 1712; C = V'o/2, 



2 ^ 



'2 - 3/32 tan 



/2^:3:5^1n^ 



(3.9) 



Let us now consider the structure of the solution in the intermediate region. We will start 
with the qualitative discussion of the non-confined solution of the type I. In this case the 
potential Vo{C,) in (|3.3b|) without pion contribution behaves as shown in Fig. |l^ (solid line) 
and the solution mo(0 corresponds to the dashed curve in the same figure, where A is the 
QCD scale at which /J^ = 2/3 and rric is the quark mass. 




Figure 12: 

14 



The pion-induced potential 



K(0 



decreases at large and small Q^ values: 



-7^/5 COS - 



(3.10) 



K 



f 



^m 



at Q 



2 ^ 2 



K 



PQ' 



at 



g2>A2 



If / ~ m ~ z/ < A, as it is natural to expect from the expression for VJr, the total potential V 
becomes essentially different from Vo{^), especially for small /5(0) values. It will correspond 
roughly to the second solid curve in Fig. |1^. The solution m(^) of the equation ( |3.3b[ ) also 
changes substantially. Its exact form depends, of course, on the values of /, u and on the 
behaviour of /3(0). At sufficiently small (3{0) it can even have a shape m(^) corresponding to 
the second dashed line in Fig. |12|. From the point of view of the structure of the solution in 
Euclidean space, these changes do not matter too much (at least at the first sight). But, as 
we will see, from the point of view of the analytic continuation the situation is different. 

Let us consider now the behaviour of the potential V{C,) and the solution u{C,) corre- 
sponding to the "confined" trajectory C. In this case Vo{C,) and Uo{C,) behave as shown in 

Fig.ra. 




^2 Q2 Qo Qi A 



Figure 13: 



The pion potential Kr(0 iii this case is 



-/3« 



m2 



2/3 



l + cos(2v/2-3/32 1n^ 



at Q > \ 
at Q < 1712 



(3.11) 



This potential is strongly localised near the value Q = Qo where the trajectory C reaches its 
highest point (point ipo in Fig- 0)' ^^ ^^^^ highest point is in the region close to zero, the 
minimum of K-(0 equals —QlP{Qo)/ P- As a result, the potential V{^) has a behaviour 
shown in Fig. |13|. The solution m(^) looks qualitatively the same. At some point its derivative 
u{^) changes sign, which leads to the transition from acceleration to damping in the equation 
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( p. 3a ). As a result, the trajectory returns to the mininium of the well. Obviously, there is 
always a solution for m(^) with such a behaviour. Because of this, a solution of the type 
C for '0((^) also exists and depends on two parameters vi and Qq. The problem with this 
solution is not its existence in Euclidean space but its behaviour at complex q^. 

In order to understand the problem of the analytic continuation, we have to remember 
that the singularity of G~^ corresponds to the point where 

u{q) = 

and find out how the position of this zero depends on the structure of u{C,) at negative q"^. 

Let us start with the simple quasi-classical solution for u: 



« = Mocoshc(o , m= [ VvW)dC 



(3.12) 



where ^o is chosen so that u{^o) = 0, and consider the analytic continuation of ( |3.12| ) along 
the circle in the complex plane of g = Qoe^^^~^\ X ^ 0; ^o = InQo (see Fig. |1^). 





Imq = Q 


Q2 


- - - , ^ 


Qo 






■^ \ 




N ^ 


Qi 






X ^ \ V 




^ X^ \ \ 




xy ' 




/ 1 1 1 J 







Reg 



In this continuation we have 



C(0 



X 



Figure 14: 



i-n 



-M \^jV yiQ + — ~ix']dx =Ci-K2, 



u = cosh ^1 cos C2 — "^ sinh ^1 sin ^2 



(3.13) 



u will have a zero at the point C2 = f ? Ci = 0. The condition ^1 = means that the 
singularity appears near the circle where u(^o) = 0. The condition C2 = f gives 



C2 = xVv{^, X 



IT 



^VvW)' 



(3.14) 



Hence, we will have a singularity in the upper half-plane if V{^q) > 1, and in the lower 
half-plane, on the unphysical sheet, ii V{^o) < 1. 



16 



From (3^) it is obvious that if the pion potential is not present, then V{^) > 1 and there 
will be a singularity on the physical sheet if the solution has a minimum (the acceleration 
transforms into damping). But with the pion potential included, ^(^o) can be less than unity, 
and the singularity will be on the unphysical sheet if we choose a solution with a minimum 
close to the minimum of the potential. At the same time, considering the behaviour of the 
solution for the non-confined trajectory (Fig. |l^), we see the opposite situation: without the 
pion potential Mo(0 has no zeros and thus there is no reason for singularities. However, if in 
the presence of the pion potential at a small /?(0) the solution m(^) shown in Fig. |l^ exists 
and has a minimum in the region where V{^) > 1, we will definitely have a singularity in the 
upper half plane. It is important to stress that for the non-confined trajectory the position 
of the minimum is defined by the form of the potential and by the condition according to 
which the solution has to approach zero in the limit Q — > 0. This means that at least for 
g{0) close to unity the non-confined solution will be unstable. 

We shall use this qualitative considerations as a hint how to obtain constructively the 
stable solution for the confined trajectory. Let us consider the trajectory C for which ipo, 
defined as the point where 'u(^o) = 0, is very small, ip <^ 1, and so is V'l where ip = 0. In 
this case the potential V{^) in ( |3.3b| ) will be determined mainly by the pion potential. In 
addition, it follows from ( p.3a| ) that ip{C,) will be a slowly changing function of C, near ,^o 



where ip = and u = 0. In this region the equation ( |3.3b| ) can be written as 

u — u'^u + x'^fx^u = (3.15) 

where u'^ = 1 + (3"^ (3 sin^ ^ + ^ j , yU^ = l3~^ cos^ ^ and x = j(3 with slowly varying u and 
yU, z/^ ^ 1 and yU^ ^ (3~^. 

For constant /i and z/ values, ( p.l5| ) is the Bessel equation in x ~ In^. Its solution can 
be given in the form 

u = CiY^{fix) + C2 J,.(/ix). (3.16) 

We choose the coefficients ci, C2 in such a way that 'u(^o) = at a point xq = Qof /Pi^o) 
inside the region between Qi and Q2 where z/ is close to unity. 

Outside this region u{C,) 7^ 0. Indeed, in the x — > limit only the singular Bessel function 
is important, 

u = ciYy{jix) ~ {fi'X)~'^ . (3-17) 



For large x, as it can be seen from (|3.11|) , u also has a power behaviour ( p.l?| ) (only with 
different fi{Q) and i^{Q))- 

Thus we conclude that in the dangerous region where u{^) can be zero we have an explicit 
expression for 

m(0 = CiFi(x)+C2Ji(x), (3.18) 



which allows us to carry out the analytic continuation along the strip shown in Fig. [T^ and 
to see where the zeros of u{^) are. 
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The most interesting case is when the singularity appears near the real axis in the q 
plane: X = f + "^^ 

Ji(e-^f x') = -ih{x') , 

2 (3.19) 

Y^{e-'^x') = --iK^{x') - h{x') . 

n 

The equation for the position of the singularity is 



—z 



-K,{x')-^h{x') 



-/i(a:') = 0, 7 = -^ = #4. (3.20) 

Cl Ji(Xo) 



It is clear from this equation that there is no singularity on the real axis. Near the real axis, 

x' = xe~^^ ^ X — iSx, we have 

-Ki{x) - 7/1 (x) - /i(x) -6 = 0, 
2"; . X . (3.21) 

-K^{x) - 7/1(0;) j ■ 5 + /i(x) = 0, / = xd^f{x) . 

If X is not small, there is no reason to expect a singularity near the real axis since in this 
case the potential in ( p.l5| ) is much less than unity. For small x values we have 

-^i(x) ^ -fl + ^ln-); /i(x)^f. (3.22) 

TT TXX \ I C 1 



The 7 parameter for a;o <^ 1 becomes 

n(a:o) 4 



7 = ^^^ = — ^ > 1 , (7 > 1 for xo < 3) . (3.23) 



Ji(xo) 



vrx, 







Using (|3.22| ) it is easy to derive the solution of (|3.21 



(5 ~ — , x^ -— ^ xl. (3.24) 

27 7r7 ° ^ ^ 

Hence, the singularity is on the unphysical sheet (5 > 0). It is positioned on the dangerous 
circle a; = xq and near the real axis when xq is small. 

Performing the analytic continuation we supposed that ^ is changing slowly as the func- 
tion of q. This hypothesis is correct, if the singularity is far away from the real axis (5 > 1). 
If, however, 5 <C 1, we have to be careful: near the singularity the trajectory 0(^) changes 
fast, jumping to a large value ~ In 1/5. 

Fortunately, this jump does not influence essentially the behaviour of m(^) or the position 
of the singularity. At first sight this jump is important because near the singularity we have 

/i q 



and in the equation for u two terms u and ucjP' /A are of the same order. Let us, however, 
consider the equation for z = v^. We can write this equation in the form 



z-2[l- 3/3'sinh^ ^ + ^/3cosh' ■^]z--\^- (3'(f)' ]z = 0. (3.26) 



2 ' ip^'''^'f)'-l{7^ 

It is easy to show that the last term in this equation has a finite hmit E when q'^ -^ /x^. 
Hence, the equation can be rewritten as 

z-2(l-3P^ sinh^ ^ + ^P cosh^ f) z + E = . (3.27) 

Its solution near the singularity is 

z = {q'-^') + c{q'-f,'f~^ ■ e^^{i?-q'r''^ (3.28) 

and the correction to ;z = m^ from the jump in the trajectory turns out to be small if /9 > |. 



4 Analytic properties of the Green's functions 
of confined quarks 

In the previous sections we have shown that the equations for Green's functions are not 
sufficient to define the theory. It is easy to show that without confinement the requirement 
of usual analytic properties for the Green's functions is enough for defining the solution for 
the Green's functions. But even in this case it is not clear a priori, what kind of bound 
states have to be introduced and how these states would infiuence the equations for Green's 
functions. In particular, we have to see how Goldstone type states which necessarily have to 
exist in some cases will affect the equations. This question will be discussed later. Having it 
in mind, we first have to learn what kind of analytic properties are necessary for the Green's 
functions of quarks and gluons if there is confinement. 

Let us start with the usual definition (given by Dyson) for the Green's function of the 
quark: 



Gd{x-x') = {T^{x)^{x')) 

= i9{xo - x'q){^{x)^{x')) - i»{x'q - Xo){^{x')^{x)). 

The quantity on the right-hand side can be written in the form 



(4.i; 



mxMx')) = I ^/(p)e-'^(--^') , (4.2) 



where 



/(p)=7or(p)7o = /(p), (4.3) 
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r{p) = J2mmn){n\^{0mS{p-pn), (4.4) 

n 

and 

mx'Mx)) = J (0i/(p)e-^^(^-^') ; (4.5) 

7{p) = Kp) , (4.6) 

hp) = Y.mmn){n\^{m)5{p+pn)- (4.7) 

n 

Making use of charge conjugation, we can see that 

hp) = C-'f{~p)C . (4.8) 

The usual hypothesis which leads to simple analytic properties for G{p) is 

pI>0, p„o > , (4.9) 



for the energy- momentum of the intermediate states pn in (|4.2|) and ([4.5|) . If this were true 
the quarks would be observable. 

If there are no stationary states with quark quantum numbers, the functions f{p) and 
f{p) are still related by (^78| ) and they still satisfy ( |4.3| ) and ( [4.6[ ) but no other simple 
properties. It can be said only that the Dyson T-product Go can be written in momentum 
space as the sum of two analytic functions Gi, G^- 

Gd = Gi-Gs, (4.10a) 

^1(^0,^1 = - r ^^^^—f{Po,P), (4.10b) 

G2iPo,p) = - r , '^^'\ . f{p'o,P) (4.10c) 

Nevertheless, we will see that the requirements of causality and unitarity impose severe 
restrictions on the functions f{p) and /(p). 

Let us regard first the causality. In order to exploit it, consider the retarded Green's 
function 

Gnix-x')=^ixo-x'o){-^ix),^{x')} . (4.11) 

In momentum space Gji can be written in terms of the same functions f{p) and f{p) 

Gr{po,P) = Gi{po,p) + G;{po,p} , (4.12) 

where 

G;{po,p) = - / . ^° . f{p'o,P) ■ (4.13) 
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The Green's function ( [4.1 1|) has to be zero for (x— x')^ < 0. This condition can be formulated 
in terms of f{p) and f{p) as 



f{p) = -Kp) at p2 < 



(4.14) 



The resuhs (|4.1CI| ) and ( [4.12D -( p:.14| ) can be expressed in terms of one analytic function in 
the complex plane with cuts as shown in Fig. |T5|. 



-p\ + ie 


® 
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-p\-ie 


p\-ie 



Figure 15: 



The Dyson Green's function is defined for real po values and has an imaginary part different 
from zero everywhere. The retarded Green's function is defined in the upper half plane and, 
due to the condition ([4.14|) , has no imaginary part between —\p\ and \p\. This means that 
Gf( can be defined in the complex plane with four cuts (Fig. |T6|) provided ie is finite. 



-p\ + ie 
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p\ +ie 
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Figure 16: 

In order to find the analytic properties corresponding to Fig. ^, the cuts denoted by the 
dashed lines have been moved from the region —\p\ < Po < \p\ to the regions po > \p\ and 
Po < ~\p\- Of course, this procedure defines a new function, different from Gd, on the real 
axis; it is natural to call it the Feynman Green's function Gp. It is, I believe, possible to 
show that this function corresponds to the Fourier component of 



Gf{x,x') = {TS^{x)^{x')) 



(4.15) 



where S is the S'-matrix in interaction representation; |), (| are free vacuum states. Gf{x—x') 
is usually calculated as a series of Feynman diagrams. Generally speaking, (|4.1|) is not equal 
to(|15D. 
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Independently of introducing Gp, Fig. |T6| suggests that we can define two functions G+, 
G_. G^{p) is the contribution of the usual cuts (|p | — ie, cxo), (— cxo, — |p | + ie) 



GM 



oo 



P'o 
dp"^ 



Po 



1 



n 



-\p\+i£ 



dp'o^+{-po,p) 
P'o-Po 



(4.16) 



'2\ 



iQ-ie P P 

and G^(p) is the contribution of the new cuts 



G-iP) =- r -^a4p'] 



(4.17) 



The expressions ( [4.16| ) and ( |4.17| ) have to be understood in the usual way, as equations for 
two invariant functions. The discontinuities A+ and A_ are now complex but A_|_ + A_ is 
real. The unitarity relations which under the usual (non-confined) circumstances define the 
discontinuities on the normal cuts have to define also the new discontinuities. To find them 
is equivalent to finding the functions f{p) and f{p) which are the imaginary parts of the 
Dyson Green's function and the retarded Green's function, respectively: 



Im Gd{p) 
ImGnip) 



f{p) - hp) , 
f{p) + hp) ■ 



(4.18a) 
(4.18b) 



Unitarity relations for quark and gluon 
Green's functions 



Let us try to calculate {^ {x)^ {x')) using the hypothesis that the interaction can be intro- 
duced adiabatically in the theory. This idea is not obvious at all. In the normal perturbative 
case it means that starting with a free theory and introducing the interaction adiabatically, 
we gradually reach the vacuum state corresponding to a theory which includes interaction 



(curve a in Fig. T7). 



^E 




Figure 17: 

But if the spectrum of the theory with interaction differs drastically from the spectrum of 
the free theory, this hypothesis will not work because there has to be level crossing. The 
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curve which has a chance to reach the physical vacuum (curve b in Fig. \[T\} starts not at 
a vacuum state with g = but at an excited state. In our case the adiabatic hypothesis 
is equivalent to the assumption that there exist some excited states in the free theory from 
which the physical vacuum can be reached by introducing the interaction adiabatically. 



In order to calculate (\I'(x)\E'(a;')), let us consider the Dirac equation for the quark field 



*fx 



(9-i)^-m^ = 0; 



(5.i: 



here A is the operator of the potential for the colour field. The operator \E' can be written 
in terms of its initial conditions at t ^ — oo in the form 



^(x) 



d'y [doGUx - y)^{y) + G%ix - y)do^{y) ] 






(5.2) 



where Gn is the retarded Green's function of the free Dirac equation. Diagrammatically 
( ^.2|) can be expressed in the form 




(5.3) 



y,yo ^ y,yo ^ y,yo 

Expressing A through its initial conditions we will have, for example, the diagram 

X 



^(x') 




(5.4) 



• y,yo 

instead of the second diagram in (|5.3| ); the dashed line corresponds to the solution of the free 
equation for the gluonic potential. Doing now the same for ^{x') with the advanced Green's 
function, we obtain for ^(a;)^(x') the following diagrammatic structure, where yo = y^. 






(5.5) 



■^o{yy '^o(y') *o(l/)' '*o(l/0 A^{yY Mo(2/' 
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Let us denote 



(^(x)*(x')) 

(A(x)A(x')) 



-^f- 



(5.6) 



In order to calculate them, we have to take all possible average values of A A and \E'\E'. As a 
result, for (\E'(x)^(x')) we will have the equation 



-X 



-^f- 



+ ■ 



X 



X 



(5.7) 



-^f- 



y ^ ~ y' 



and a similar graphic equation for {A{x)A{x')). The functions corresponding to the vertical 
lines are now exact retarded Green's functions 



, 'X- 



G 



R 



Gl 



G'n 



-^f- 



Dl 



+ 



(5.8) 



Let us see, what happens to ( |5.7|) when yo = t ^ —oo. In a theory without confinement the 
answer is very simple. If the retarded Green's function has a pole at p^ = m? , the diagram 



x 



y 



-^f- 



X' 



(5.9) 



has a finite limit when yo ^ —oo and equals 
dqQe 






iM.-.')_d^ 5(A;2 _ ^2) ^((^ _ f^y _ ^2) ^^^) 



where f{k) corresponds to 



-^f- 



. The product of the (5-functions gives 



6{{q-kf-m')6{k'-m') = 6{k^-m')6{ql-2koqo) ^ ^(^' ""'^ ^(go) 

Z/to 



(5.10) 



(5.11) 



and therefore (|5.10| ) does not depend on yo. For a free Green's function 



X 



-^f- 



X' 



(5.12) 



and instead of (|5.7|) , we have the usual relation for the imaginary part of the fermion Green's 
function 



fip) 
->< — 



Up) 
->< — + 



'X~. 



fiP) 



+ 



(5.13) 
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Repeating the same calculation in a theory in which the Green's functions of quarks and 
gluons have no poles but soft singularities we will get (\l/(a;)^(a;')) = in the t — > — cxd limit. 
This, of course, cannot be the correct answer — in the theory with confinement there exist 
stable hadrons in general, and those belonging to the pseudoscalar octet in particular. The 
problem can be resolved by introducing these pseudoscalar particles from the very beginning 
as elementary objects [6]. In this case we include in the Dirac equation for \1/ not only the 
colour field A but also the pseudoscalar field ip. This will add to ( p.7| ) diagrams of the form 




+ \ 




(5.14) 



where the wavy lines correspond to pseudoscalar particles; for the sake of simplicity, we shall 
call them pions. They satisfy the equality 



X 




X 



X x' 



(5.15) 



Instead of ( |5.13| ), we have 



-^ — 



f\v') 



and for the gluonic imaginary part, correspondingly. 



f{p') 






+ 



(5.16) 



-X 



'-^Xr' 



X--r 



-^ 



^-^Xr-*"^ 



+ 



(5.17) 



The equation (|5.16|) differs essentially from (|5.13| ). It contains no driving term (like /q in 
( p.l3| )); it is a non-linear but homogeneous equation. The equation ( |5.17| ) contains driving 
terms coming from fermions. 



Let us consider in detail ( |5.16|) in terms of the equation ( |4.1CI|) for f{p). We can write 
f{p) in the form 



f{p) = Hp'nuip) + f-ip)} + H-p')fEuci{p) 



(5.1^ 



where f± are the contributions of positive and negative frequencies po- 

Let us first investigate the contributions of the functions /+ and /_ to the first term 



in the right-hand side of (|5.16| ). The function f+{p) gives the usual contribution to the 
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imaginary part of the quark Green's function coming from the intermediate [positive energy] 
TT and quark (g+) states. The contribution of f-{p) is very interesting because the signs of 
the particle energies in the intermediate states are different [ir and a negative energy quark 
g_]. The expression for this contribution (without the external legs) is 

A,.,. = 9' J 0[^«-(?') + b^q'm-qomp - qf - ^i'mpo - go). (5.19) 

Here g is the coupling constant connecting the vr-meson to quarks; a_ and 6_ are two invariant 
functions for the imaginary part of the quark Green's function. 



In order to understand the structure of (|5.19| ) let us consider the integrals 

8^ 



hip, K, /i) = / ^^5(g^ - K^)6{{p - qY - /i^)?9(-go)^(po - go) 



X '&(p'^)^{{k - ^if - p2){^(-po)^(K - /i) + ^(po)^(/^ - ^^)} 
and 

Ii{p,K,ix)= / -^q5{q'^ - K'^)5{{p-qY - f?)'&{-qo)i3{po-qo) 

= ^jM + «:' - /i^)V[(/.-^)2-p2][(^ + ^)2_p2] (5.21) 

X ^{j)'^)d{{K - ^f - p2) {^(_po)^(fi; - /i) + ^(po)^(/i -/«)}. 

Hence, 

Ag_,, =/ / dK^{a^{K^)h{p,K,^) + b_{te)l2{j),K,^)}. (5.22) 

The structure of Ji and I2 implies that /_ on the right-hand side gives contributions to both 
/_ and /+ on the left-hand side. The equations for /_ and /+ can be written in the following 

way: 

/_ = Gr/\-^_.^Gr + ... (5.23) 

/+ = GrA+_,.Gr + GrI\+^^^Gr + --- , (5.24) 

where A-__^(p) and A+ _^(p), 

d^[a^{K^)h{p,K,ii) + h^{K^)h{p,K,ii)]+... , (5.25) 



2 



X 



I (i/t^ [a_(/t^)/i(p, K, /i) + 6_(fi;^)/2(p, K, /i)] + . . . , (5.26) 

Jo 



are the contributions corresponding to the quark masses k larger {k > /i) and smaller [k < /i) 
than the vr-meson mass /x, respectively. Higher order terms have the same structure: large 
quark masses of /_ contribute to /_, small quark masses contribute to /+. 
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The result of this analysis is the following. The imaginary part of f-ijP') which equals 
zero in a non-confined theory has to satisfy a non-linear homogeneous equation. Due to 
asymptotic freedom, it must decrease at large p^. In linear approximation the equation for 
/_ with zero boundary conditions at infinity is a typical bound state equation defining the 
quark distribution in the vacuum. 

The equation for /+ is very different. The solution of the equation for /_ gives an 
inhomogeneous term for the equation for /+. The latter becomes similar to the equation for 
/+ in a non-confined theory ( p.l3| ). Now, instead of the usual /o(p) = ^(jP' —w?){p + m), we 
have 

Up) = GnAlJp)Gn. (5.27) 



Note that, according to ( ^.26|) , /o(p) is proportional to '&{fi — a/p^)- Iterating, we have to 
find /+ which corresponds to perturbation theory at large p^. 

The equations for /+ and /_ contain Gr. If we were able to find /+, /_ for a given Gr, 
we would then obtain an equation for Gr in the form of the dispersion relation 

Gnip) = - f -^fip'o.p) --I -^fip'o.P) ■ (5.28) 



This Gr{p) has no poles since there are no (5-functions in /_(. and /_. The equations ( |5.16| ), 
( ^.17] ) can be considered as generalised unitarity relations. 

The structure of the equations for gluonic fg+, fg_ is different only in the sense that fg_ 
contains an inhomogeneous term corresponding to intermediate quark states. 

All these equations are in excellent agreement with the idea that when the negative 
energy states are not completely occupied, both positive and negative energy states become 
unstable. The decay modes of these states are explicitly written in the equations ( |5.16| ), 

(EI3)- 

Up to now we ignored the contribution fEuci{p)'&{—p'^)- Let us include it in the right- 
hand sides of ( [5.16| ), ( |5.17| ). The equation for /_ will remain the same homogeneous equation 



(since it is impossible to have a time-like vector p with a negative po as a sum of a space-like 
quark momentum p' and a time-like pion momentum k with a positive ko). The equation 
for /+ (|5.24|) will acquire an additional term: 

/+ = Gn A+_,. Gr + Gn A^^^ Gr + Gr A+^^^^,, Gr . (5.29) 

The equation for JeucI can be written in the form 

fEucl{p) = Gr Af_";' Gr + Gr Af^^J Gr , (5.30) 

containing a driving term. The iterative solution of (|5.3CI|) has the following diagrammatic 
structure: 







/^^ = ^/^\^ + ^X>>^ + ... . (5.31) 



vr 
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Inserting this solution into the last term of ( |5.29 ) we will have an equation similar to ( ^.24 ) 



with the only difference that A j"_ ^ in ( |5.24| ) acquires an additional contribution from multi- 
pion states. This will not change any of the conclusions we have reached, because it is only 
/+ and /_ which matter for the calculation of the retarded Green's function. 

The structure of fsuci written in ( ^.31j ) is, nevertheless, interesting in two respects. 
Firstly, it shows how the analytic features corresponding to Fig. |15| (with functions hav- 
ing discontinuities at positive and negative p^ values) can be reduced to analytic features 
corresponding to Fig. 110 where the discontinuity differs from zero only at positive p^ values. 



Secondly, the existence of the integrals in (|5.31| ) imposes a restriction on f~.{p)- 
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